Effects of non-magnetic impurities on spin-fluctuations induced superconductivity 
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We study the effects of non-magnetic impurities on the phase diagram of a system of interacting 
electrons with a flat Fermi surface. The one-loop Wilsonian renormalization group flow of the angle 
dependent diffusion function D(6i, 62,63) and interaction C/(#i,#2,#3) determines the critical tem- 
perature and the nature of the low temperature state. As the imperfect nesting increases the critical 
temperature decreases and the low temperature phase changes from the spin-density wave (SDW) 
to the d-wave superconductivity (dSC) and finally, for bad nesting, to the random antiferromagnetic 
state (RAF). Both SDW and dSC phases are affected by disorder. The pair breaking depends on 
the imperfect nesting and is the most efficient when the critical temperature for superconductivity 
is maximal. 
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Interesting and new physics appears when density wave 
(DW) and superconducting (SC) correlations interfere 
with each other and with the disorder. One finds realiza- 
tions of such interferences in correlated electronic systems 
whose Fermi surface (FS) has at least imperfect nest- 
ing properties. The Bechgaard salts (TMTSF) 2 X and 
the Fabre salts (TMTTF) 2 X have two open quasi-one- 
dimensional Fermi sheets responsible for DW correlations 
and, very probably, for superconductivity. Q In these 
compounds the disorder strength can be tuned by the 
cooling speed or by x-ray irradiation, while the imperfect 
nesting is controlled either by chemical structure or by 
pressure. Another example of correlated system with im- 
portant disorder effects are the two-dimensional (2D) or- 
ganic compounds from the BEDT-TTF (ET) family. 0,0 
All ETs have a metal-insulator and/or a DW transition 
in the phase diagram. They have two FSs with several 
more or less perfect inter- and intra-band nesting wave 
vectors. Their superconductivity is induced by the DW 
fluctuations. In the high-T c superconductors the pseu- 
dogap and the superconducting states are both related 
to the antiferromagnetic (commensurate SDW) insulat- 
ing phase. The effects of disorder on the high-T c su- 
perconductivity has been widely studied experimentally 
and theoretically. 0, It has been shown by a number of 
methods, and in particular by the angle resolved renor- 
malization group (RG) that this interplay between the 
DW, the pseudogap and the SC exists in models with 
at least flat segments of the FS or with nested van Hove 
singularities. [20,111 However, no attempt has been made 
to implement both disorder and correlations in the N- 
patch RG, which would provide a perturbative, but non- 
biased and controlled way for constructing the phase di- 
agram. 

In this letter we report a perturbative angle-resolved 
(N-patch) RG theory for the model consisting of two flat 



FS segments in regard. The imperfect nesting parame- 
ter is treated approximately, by introducing a cut-off e in 
the Peierls channel. For the pure case this problem was 
solved in the Parquet formulation by Zheleznyak et al Q 
and later closely reexamined in the field-theoretical RG 
language. While previous theories in more than one di- 
mension (ID) were based on the assumption of one domi- 
nant order parameter . pi fiol Hll | we take into account in a 
non-prejudiced way both Cooper and Peierls correlations 
together with the diffusion processes, generalizing to 2D 
the ID theory of Giamarchi and Schulz.^l 

We first consider a non-interacting isotropic disordered 
electronic system with bare dispersion £(fc). We con- 
strain the momenta to lie within the shell of ±A around 
the Fermi surface and follow the Wilsonian one-loop flow 
starting from A . At the scale I = ln(A /A) the propa- 
gator has the form 



Q l - l {k>uj)=iuf l {u)-fi{k)> 



(1) 



where // is a flowing junction of u>. Considering only the 
self-energy correction due to the disorder, see Fig. QJa) 
and discarding for simplicity one-loop graphs renormal- 
izing the diffusion D (Figs.^b) and (c)), we find 
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(«fA) 2 + Ki(^)F 



with fi = o(w) 



(2) 

where r = vf/2D is the elastic scattering mean free time. 
The Born approximation /( B ) = 1 + 1/(2|cc?|t) is recov- 
ered in the limit A — > 0, when setting fi(u>) to its bare 
value fi = 1, in the RHS of Eq. Q. In Fig. we see 
how the ^-dependence of / is constructed successively 
by mode elimination. For any finite A, the one-particle 
propagator is a regular function of frequency. The charac- 
teristic momentum scale for disorder is Ad = 2D/(irvp). 
In the Born approximation, and for A/ Ad -C 1, the slope 
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FIG. 1: Self-energy (a), particle-hole (b) and particle-particle 
(c) graphs taken into account. T symbolises both interactions 
U and disorder D. Barred propagators are on-shell at A. 



of w// B) (w) at the origin is // B) (0) ~ A d /A. The ex- 
act solution of Eq. Q is much steeper since its slope is 
/*(0) ~ exp(A d /A). 

We now consider the disordered and interacting system 
whose FS consists of two flat segments in regard. Wc 
introduce left (L) and right (R) angle-resolved movers, 
and decompose the momentum k into one component 
k orthogonal to the FS and the other k\\ parallel to 



the whole frequency dependence of the propagator into 
account. Here we consider the problem from the weakly 
coupled Fermi liquid side. We want to study how the 
model is driven away from its Fermi liquid fixed point by 
the combined effects of weak interactions and disorder. 
For this purpose the zero-order scaling arguments 0] 
allow us to keep only terms linear in energ y in the non- 
interacting part of the replicated (see e. g. |llj]) effective 
action, which then reads 



So = ~ / -7T- 



dX 



iu>zr 1 {e) - v F k ip(x)ip(x). (3) 



X = (k,uj,8,a,a) is a multiple index. The integral over 
X is a shorthand notation for integrals over k, lo and k\\ 
and discrete summations over the replica index a, spin 
index a and the left-right index. In Eq. J3J the quantity 
the FS. The angular position on the FS is denoted by (9) = fi(uj = 0,8). It should not be confused with 



6 = (R/L,/c||). We will assume the Fermi velocities 
do not depend on 8. In the one-loop approximation, 
the one-particle propagator at scale I is Qf 1 (k,uj,8) = 
iLofi(uj,8) — v-pk. A complete description of the disor- 
dered phases would require a fully functional RG, taking 



the one-particle weight, discussed in In fact, the 

effective action @ can be interpreted as the marginal 
part of the diffusion pole as "seen" by the electrons at 
the RG scale I. The interaction and the disorder parts of 
the replicated effective action read 



S-nt = 27TOF 



dX l ... dX 4 

(2^) 4 



[27n5(4 + 3 - 2 - l)]6 aua J a2ta3 6 aiia3 Ui(di,8 2 ,8 3 )ip4ip3ip2'ipi, 



Sdis = 2wv F - 



v F A f A dX x ... dX 4 



(4) 



(2tt)* 



[2tt8(A + 3 - 2 - l)]6 ai , ai S a2 , a3 [2irS{ui - u A )]D,(6i, 8 2 , 8 z )M^i- (5) 
I 



We have denoted ipi = ^(X l ). The 5(4 + 3-2-1) 
function ensures momentum and frequency conservation. 
U couples electrons in the same replica only. Wc keep 
only the angle dependence of U, and drop momentum 
and frequency dependence since we restrict ourselves to 
the marginal part of the interaction. We do not take 
Umklapp processes into account because we suppose that 
the system is away from half-filling. The disorder part 
of the effective action is the gen eralisation of the one 
used by Giamarchi and Schulz.|12j Two angle dependent 
diffusion functions (forward) and (backward) are 
gathered in a single diffusion function 0(81,82,83). The 
disorder couples electrons from two different replicas and 
in the same replica. The <5(a->i — W4) function ensures the 
elastic character of the diffusion function. We keep only 
the relevant part of the diffusion function and neglect 
all marginal contributions, which are linear in momenta 
perpendicular to the FS and in frequencies. 

We derive one-loop RG flow equations by successive in- 
tegration of fast modes. The Feynman graphs generated 
by infinitesimal mode elimination are shown on Fig. ^ 



Schematically, the angle resolved flow equations read 
diUi{8i,8 2 ,8 3 ) = p u [Ui,Di](8 1 ,8 2 ,8 3 ), (6) 
8^(8^82,83) = (D l +l3 D [U l ,D l ])(8 1 ,8 2 ,8 3 ). 



with Ui(9 u e 2 ,e 3 ) = Vgggag i Ui(8 u e 2 ,6 s ) and 
A(0i,02,0 3 ) = eMl)VZiZ 2 Z 3 Z 4 0^,82,83), where 
Zi = Zi(8i). Note that the equation for the disorder 
contains a linear term in D, explicitly showing the rel- 
evance of the disorder. The /3-functions are homoge- 
neous quadratic functions in U and D. The disorder 
is not renormalized by a term quadratic in U and the 
interactions are not renormalized by a term quadratic in 
D. Indeed, starting from a pure situation, interactions 
alone cannot induce disorder and starting from a non- 
interacting disordered situation, elastic couplings (the 
disorder) cannot induce inelastic couplings (the interac- 
tions). We checked our RG equation in two limits. In 
the ID limit, neglecting the contributions quadratic in 
the disorder, we recover the (weak coupling) equations 
of Giamarchi and Schulz. 121 The effect of disorder on a 
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FIG. 2: Progressive building of the ^-dependence of fi(w) as 
modes are eliminated, (a) is the Born approximation / ; (lo) 
and (b) is the exact solution of Eq. J5J 
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FIG. 3: Effect of disorder on the phase diagram. Inset shows 
the critical pair breaking parameter. The dip at e = T c q 
corresponds to the maximal SC critical temperature for the 
pure case. 



conventional s-wave superconductor is obtained from our 
equations when dropping the particle-hole contributions 
and using initial conditions of the attractive Hubbard 
type. In this case we recover Anderson's theorem. |l5| 

After discretization of angular variables, we have in- 
tegrated numerically the equations ©. The robustness 
of our results upon increasing number of patches is sat- 
isfactory above N — 16 patches per side of the FS. The 
dependence of the critical scale on the imperfect nesting 
gives the phase diagram. For the pure case, i.e. the pair- 
breaking parameter a = 1/r = 0, we reproduce the result 
of Zheleznyak et al, shown in dotted line on Fig. [3J As 
one can conclude from the flow of dominant eigenvalues 
of the coupling function U , shown in thin lines on Fig.^J 
the part of the phase diagram with T c > e shows a tran- 
sition to the SDW state with enhanced dSC fluctuations, 
while the T c < e part corresponds to the onset of d-wave 
superconductivity. We now introduce a finite disorder 
via the pair breaking parameter a = 2Dq/vy — 0.39T c o, 
Dq being the bare diffusion constant. The SDW and dSC 
phases are both affected by disorder. For e grater than a 
disorder dependent critical imperfect nesting e* the dSC 
phase is completely suppressed and a new low tempera- 
ture state appears. The signatures of all three phases are 
visible on Fig.0] At perfect and good nesting, i.e. as long 
as e < T c , the flow near the divergence is almost the same 
as for the pure case: the SDW correlations are dominant 
while the dSC ones are also divergent but weaker. The 
disorder remains finite at the critical scale. For moder- 
ate imperfect nesting, T c < e < e*, the dSC becomes the 
dominant instability because the particle-hole logarithms 
are cut. The disorder remains finite, but is stronger than 
in the perfect nesting case. Finally, at very bad nest- 
ing, e > e* the most diverging "coupling" is the disor- 
der. The dSC is irrelevant, while the SDW amplitude 
shows the intriguing re-appearance at the very last mo- 
ment before divergence. This divergence is induced by 
the particle-particle U-D terms in the renormalization 
of the interaction U. Its very small residue reflects the 



fact that the SDW interactions which are driven by dis- 
order are only these with total momentum equal zero. 
Consequently only a small subset of all SDW amplitudes 
diverges. Altogether the flow in Fig. 0Jc) allows us to 
characterise rather precisely the corresponding low tem- 
perature state : it is a charge-localized state with weakly 
enhanced SDW correlations. We expect that the proper- 
ties of this phase are similar to the ID random antiferro- 
magnetic (RAF) phase, discussed in 

We have calculated the critical pair breaking a* as a 
function of the imperfect nesting parameter e. The re- 
sult is shown in the inset of Fig. [3] As long as the flow 
is entirely in the parquet regime the pair breaking ef- 
ficiency increases with imperfect nesting, i.e. a*/T" =0 
decreases. At the point e = T c , where the flow starts 
to be affected by the imperfect nesting, a* /T c a=0 reaches 
its minimum. This is the point of the most efficient pair 
breaking. It corresponds to the maximal value of the 
dSC critical temperature. As e increases further the pair 
breaking gets progressively less efficient. The dependence 
of the pair breaking on the imperfect nesting is a very 
consequence of the interplay between Pcicrls, Cooper 
and localization tendencies through the angular depen- 
dence of corresponding amplitudes. In the Abrikosov- 
Gor'kov theory, ^^,^3 where particle- hole logarithms are 
neglected, and a dSC order parameter is phenomcnolog- 
ically assumed, the ratio a*/T c a=0 ~ 0.88 is independent 
of nesting. 

In conclusion, we have used the one-loop N-patch 
renormalization group to compute the non-magnetic dis- 
order effects on a system with spin-density wave fluctua- 
tions induced d-wave superconductivity. The pair break- 
ing is the most efficient when the T c for superconductivity 
is maximal. 
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FIG. 4: Flow of dominant components of coupling U and diffusion D, for (a) perfect nesting, (b) imperfect nesting (l e 
ln(Ao/e) =— 8.5) and (c) bad nesting (l e ~ 6.8). 



and discussions. 
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